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Abstract 

The recently released Planck 2013 results show that the primordial fluctuations are deviated from 
isotropy. The deviations from isotropy are robust against component separation algorithm, mask 
and frequency dependence. To incorporate the Planck's data into standard cosmological model, 
we propose an inflation of the very early universe in an anisotropic spacetime. A generalized 
Friedmann-Robertson- Walker (FRW) metric is presented in the Randers-Finsler spacetime. By 
employing the osculating Riemannian approach, we obtain the primordial power spectra of the 
scalar perturbation with direction dependence, such as the dipolar modulation. This is consistent 
with the deviations from isotropy of the universe found by the Planck satellite. 
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I. INTRODUCTION 



The cosmological principle assumes that the universe is statistically isotropic and ho- 
mogeneous at large scales. Recently, this principle has been rigorously tested by the CMB 
anisotropy observations from the Planck satellite l|. The Planck satellite [2I has found 
deviations from isotropy (around 3a). Several anomalies of the CMB anisotropy have been 
confirmed, such as the quadrupole-octopole alignment js], the hemispherical asymmetry 
0, Isl, the parity asymmetry j6-12|, and etc.. Besides the Planck's results, the Wilkinson 
Microwave Anisotropy Probe (WMAP) 13|] has previously found evidence for the devia- 
tions from isotropy and the anomalies. Via dealing with the WMAP dataset, the level of 
statistical anisotropy was constrained to be (7 ~ 0.15 in the primordial power spectra of the 
form 'P(k) = V{k){l + g(k- n)^) where n denotes a privileged direction in the space 14j-|l6|. 
These anomalies might stem from a non-standard spinor or certain primordial vector 

fields \mm\. 

To incorporate the observed deviations from isotropy in the Planck data into the standard 
cosmological model, we should try to investigate properties of inflation at an anisotropic 



spacetime. The natural framework of anisotropic spacetime is the so-called 



Finsler geometry gets rid of the quadratic constraint on the metric 



insler geometry 



29|-|3l|. The Finsler 



spacetime admits certain privileged axes and permits less symmetries than the Riemann 



one 



32|-l34j| . Thus, the Finsler geometry is a reasonable candidate to reveal the deviations 



from isotropy of the spacetime. For instance, the 



35 



insler spacetime could account for the 



43| . The cosmic acceleration could be 



Lorentz violation as well as the CPT violation 
explained by the anisotropic Friedmann equation in Finsler cosmology 
scale bulk flow 



45 



M|. The large 



46| could be revealed by a Finslerian Zermelo navigation model |47 |. 



The Randers structure could account for the privileged direction 



accelerating expansion in the Hubble diagram 



48l-l50l| of the maximum 



5l| . A spatially anisotropic Finslerian model 



could account for the mass discrepancy problem of the Bullet cluster [53|, and etc.. 

The primordial vector fields may be related with the Randers spacetime. The Randers 
structure 5^1 involves an extra 1-form which is related to a vector field. This vector field may 
influence the very early evolution of the universe if it existed in the universe. In this paper, 
we propose a generalized FRW metric in the Randers spacetime. The Randers structure 
comprises the FRW part and the weak vector field. The vector field singles out a privileged 



axis in the very early universe. We study Einstein's gravitational field equations with the 
osculating Riemannian approach. The time-time component of Einstein's field equations 
is resolved to obtain an inflationary phase of the very early universe. The anisotropic 
modifications are studied for the inflationary universe. By studying the equation of motion 
for the inflaton field, we obtain the primordial power spectra of the scalar perturbation 
with statistical anisotropy. The predicted statistical anisotropy may account for the Planck 
observed deviations of isotropy and certain anomalies in the CMB anisotropy. 

The rest of the paper is arranged as follows: In section [ITl we first present a brief introduc- 
tion to Finsler geometry. In section [XTTt the osculating Riemannian approach is used to study 
a generalized FRW metric in the Randers spacetime. The geometric objects are obtained. 
In section llVt we study the time-time component of Einstein's gravitational field equations 
and resolve it to obtain an inflationary solution. Meanwhile, the anisotropic modifications 
are acquired for the inflationary phase of the universe. The primordial power spectra with 
statistical anisotropy is got through the equation of motion for the inflaton field in section 
IVl Conclusions and discussions are listed in section |VT1 



II. THE ANISOTROPIC SPACETIME 

In this paper, we suggest the Randers spacetime as a suitable background of inflation at 
the very early stage of the universe. The Randers spacetime is a class of Finsler spacetimes. 
Thus, we first present a brief introduction to the Finsler geometry in this section. For more 



detailed discussions on Finsler geometry, see the references, for instance |29l431|. 

Different from Riemann geometry, Finsler geometry is defined on the tangent bundle. 
Finsler space stems from the arc-length integrals of the form 

F{x,y)dT, (1) 

where x and y = dx/dr denote the location and the velocity, respectively. The proper length 
is given by r (or s). The integrand F{x,y) is called the Finsler structure. It is a smooth 
and positive function on the tangent bundle. It is positively homogeneous of degree one, 
i.e., F{x, Xy) = XF{x,y) for A > 0. The Finsler metric is defined by the Hessian, 

d d ri 



Together with its inverse g^'^, it is used to lower and raise the spatial indices of tensors. For 
the Finsler spacetime, the spatial indices run from 1 to 3 and the temporal index runs in 
this paper. 

The Finsler metric is dependent on the positions x as well as the directions given by the 
fibre y. This is different from the Riemann metric which depends on x only. Thus, there is 
other significant object, i.e., the Cartan (torsion) tensor in Finsler geometry. The Cartan 
tensor is defined as js^ 

~ 2dy- • ^ ^ 

It completely characterizes the deviation of the Finsler space from the Riemann space. The 
Finsler space becomes Riemannian if and only if C^^„ = 0. This result reveals that Finsler 
geometry is a natural generalization of Riemann geometry. 
The Randers space 



54| is a kind of Finsler spaces. It could be viewed as a Riemann 
space influenced by a primordial vector field (for instance, the electromagnetic field). The 
Randers structure is given by 

F{x,y) = a{x,y) + I3{x,y) , (4) 

where a{x,y) = ^Ja^y{x)y>^y'' denotes a Riemann structure and l3{x,y) = bi^i_{x)y'^ is a 1- 
form. However, the vector field appears as a part of the space structure in the Randers 
space. Furthermore, it induces the anisotropic properties of the Randers space. The tildes 
here denote raising (or lowering) the indices with the Riemann metric a^^, (or a'^'^). To reveal 
the anisotropy of the Randers space, one notices that the Randers structure is not absolutely 
homogeneous of degree one. Otherwise, the Randers structure reduces back to the Riemann 
one. The reason is that a would keep invariant while /3 changes its sign under the rescale 



y — > ~y. In addition, the anisotropy of the Fins 



number of solutions of the Killing equations 



32|-|34| 



er spaces could also be revealed by the 



In the Randers space, the Randers metric is defined as [30 1 



dlJ.ll ( ^^^U ) ~l~ ^fM^U ) (5) 



a 



where = + and = a^^y^ /a. Its inverse is given by 



where = b'^bfj,. The Cartan tensor of Randers space is given as 30 1 



where Sf^^u^^ refers the summation over the cychc permutation of indices. We note that each 
term in C^j^a- is proportional to the components of b. Thus, the Cartan tensor could be a 
first order object in the case of <^ 1. 



III. THE OSCULATING RIEMANNIAN APPROACH 

In this section, we consider a generalized FRW metric with weak anisotropy in the Ran- 
ders spacetime. Following Stavrinos et a/.'s approach js?], we will employ the osculating 
Riemannian metric to approximate the Randers metric. 

The related Randers structure is given by the spatially flat FRW metric d^^, = (1, — a^(t), — a^(t), — a^(t)) 
and an extra 1-form b^dx^ (or a vector field b'^d/dx^). The temporal coordinate denotes the 
proper time measured by a comoving observer while the spatial coordinates are comoving. 
The velocity = {1, 0, 0, 0) denotes the tangent 4- velocity of the comoving observer along a 

6|l < 1 in the 



55|. Thus, we 



worldline. Here r denotes the proper time. The vector field is very weak, i.e., 
Randers spacetime. It would be expected to align with the tangent 4- velocity 
consider the timelike vector field that has only the temporal component non- vanishing, i.e., 
bfj, = (-B(z), 0, 0, 0). In addition, the component B{z) is set to be dependent on the third 
spatial coordinate z only. This proposition would induces the anisotropy of the primordial 
power spectra along the z-axis. Therefore, we will study the generalized FRW metric (i.e., 
the Randers structure) as 



dr = ^dt^ - o?{t) {dx^ + dy^ + dz^) + B{z)dt , (8) 

where the spatially flat FRW part is given by the first term on the right hand side. The 
anisotropy of the obtained Randers spacetime is characterized by the second term on the 
right hand side of the above equation. For the weak vector field, the Randers metric could 
be viewed as a slight modification to the FRW metric. 

Usually, it is difficult to discuss the gravitational issues in the completely Finsler geometric 
framework. However, the Finsler metric could be approximately related to the Riemann 



metric in certain cases 



29 



561 ]. Actually, it is convenient to study the gravity with the 



osculating Riemannian method 29(]. The Finsler structure corresponds to the osculating 
Riemannian metric as 

9,iu{x) = g^u{x,y{x)) , (9) 

where the velocity y{x) is viewed as a function of the position x. Then the osculating 
Riemannian structure becomes ds"^ = (1 + B^z))"^ dt^ — a'^(t) (1 + B{z)) (dx"^ + dy"^ + dz^) 
for the above generalized FRW metric ([8]) in the Randers spacetime. Correspondingly, the 



Christoffel symbols are defined as |57| 



where 7^<^(a;, = ^ - ^ + ^). From the equation ([7]), we may note that 
the Cartan tensor is proportional to the components of the vector h which is small. In 
addition, we have chosen the comoving coordinates in the Randers spacetime. Thus, the 
terms dependent on the Cartan tensor could be dropped in the equation ( ITOl) . In the Randers 
spacetime, therefore, we approximately obtain the Christoffel symbols as 

rtv(x)^7^;.(a;,i/(x)) . (11) 

Here r(^^(x) denote the coefficients of the osculating Christoffel connection. 

Corresponding to the osculating Christoffel connection, the components of the curvature 
tensor are given by 

p A» =y^^ _ F^^ _|_ F*^ F^' _ F'* F^' fi o'l 

per vo^p yp^o V(7 Kp up Ka ' \^ J 

The Ricci tensor and the scalar curvature are, respectively, given as 

Ric^i, = , (13) 

S = g^''ix)Rtc^, . (14) 

The osculating Riemannian metric evolves following the conventional Einstein's gravitational 
field equations, 

RiCf,:, - ^gpuS = SnGT^^ , (15) 

where T^,^ = d^cpdyCp — g^y {\ g"^ da4>dj3(f) — V{(f))) denotes the energy-momentum tensor 
of the cosmic fluid. However, the anisotropic properties of the Finsler spacetime have been 
approximately comprised into the osculating Riemannian metric. Thus, the obtained Ein- 
stein's field equations are weakly anisotropic in the Randers spacetime considered in this 
paper. 



IV. THE INFLATIONARY PHASE WITH WEAK ANISOTROPY 



In this section, we study Einstein's gravitational field equations in the Randers spacetime. 
The inflationary solution is obtained for the very early universe at the zero-order assumption. 
Further, we take account the slight modifications to the inflationary phase such that the very 
early spacetime acquires the weak anisotropy. 

First, we calculate the time-time component of the Einstein tensor to study Einstein's field 
equations in the Randers spacetime. The Einstein tensor is defined as Gfj_y = Ric^^ — \g^uS. 
We obtain its time-time component as 

where the dots denote the derivative of a with respective to t while the primes denote 



the derivative of B with respective to z. For an inflationary phase |58N62| . the very early 
universe undergoes a process of exponential expansion, i.e., a ~ e^* where H is the Hubble 
horizon. Thus, the second term on the right hand side of the equation (fT6l) would decrease 
exponentially. It could be discarded at the zero-order approximation. In this way, Einstein's 
field equation would reduce back to the conventional one. To account for the spacetime 
anisotropy, we could consider the modifications from the second term on the right hand side 
of ([16]). 



At the zero-order approximation, the time-time component of the Einstein tensor becomes 

= 3 {a/af . (17) 
On the other hand, the energy-momentum tensor of the cosmic fluid is characterized by the 



one of the perfect fluid. At the zero order, its time-time component is given by 63 1 



Tr=p(°) = i(^)' + V^(°H0(°)), (18) 

where 0'^°)(t) is the zero-order part of the mostly homogeneous inflaton field x). Consider 
the slow-roll condition that the inflaton field slowly rolls down its potential, i.e., (^^^-^ ^ 
y(o)(^0(o)). Thus, the energy density of the inflaton field almost becomes a constant. 
The time-time component of Einstein's equations becomes 63 1 



f^y = ^l^W(0(°)) = const. , (19) 



which determines the evolution of scale factor a. This equation has an exponential solution 
a{t) ~ c^* where H denotes the Hubble horizon which is almost constant. The exponential 
expansion of the very early universe corresponds to the inflationary phase of the universe . 



63| 



Note that the above analysis is similar to that in the standard inflationary model, see 
for details. 

To account for the anisotropic effects, we substitute the exponential solution a(t) ~ e^* 
back into Einstein's field equation (fT5|) and obtain the anisotropic modification B{z) of the 
very early universe. In this case, the time-time component of the Einstein tensor is given by 
the equation (1TB|) the energy density of the infiaton field is given by Tqo = | j + (1 + 
Here the potential of the infiaton field is set to be = (1 + Bf, 

where the correction (1 + 5)^ in Too comes from the osculating Riemannian metric goo- Then 
the time-time component of Einstein's field equation (fTSll becomes 



35'2 - 45" (1 + B)=0 , (20) 



where we have used the equations (fT6|) . (fT7|) . (fT8|) and (fT9l) in calculations. We see that 
the remained term in the equation (l20l) determines the anisotropic modifications to the 
inflationary phase of the universe. The equation (!20l) has a solution as 



BW = ^tac.)*(l + |)^l. (21) 



where Ci and C2 are the integral constants. In addition, C2 > denotes a given distance 
scale which could be constrained by cosmological observations. We could set |ciC2| ~ 4 and 
\z\ ^ C2 because of the condition <^ 1. Note that B{z) would be a monotonically 

increasing function of z. To the second-order approximation, the above solution B{z) could 
be expanded asB{z) — fo + fiz+ g(^-^^^j^_^y2 z'^ +o{z'^) where fo = |ciC2| — 4 and /i =4(H-/o)/c2. 
In this way, we obtain the inflationary phase of the generalized FRW metric ([8]) with weak 
anisotropy in the Randers spacetime. 
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V. THE PRIMORDIAL POWER SPECTRA WITH THE DIRECTION DEPEN- 
DENCE 



In this section, we study the primordial power spectra 64j-l70| with weak anisotropy in 



the Randers spacetime. First, the action of the inflaton field is given by 

S = J (fx^gC = j d^x^ (lgf^'^d,(f>d,(j) - Vicj))^ (22) 

where y/—g = det{g^^) for the osculating Riemannian metric of the Randers structure ([8]). 
From the Euler-Lagrange equation ^§ - d, {y^^^ = 0, we obtain the equation 
of motion of the infiaton field as 



--^-i^VV.^-|5.-0. (.3) 



We have decomposed the infiaton field into two parts 4>(t,x) = + (50(i(:,x), where 

S(j){t, x) denotes a first-order perturbation. 

For the zero-order part ^'•^^(t), the equation of motion becomes 

0(0) + 3if0(o) + ^^^^^^ = . (24) 

For the anisotropic and inhomogeneous perturbation, the equation of motion for the fiuctu- 
ations becomes 



^2t/(0) 1 _|_ R R' 



where we have used the equation ( l24l) to eliminate Typically, the ^^^5— term is small, 
which is proportional to the slow-roll variables. Thus, it could be neglected in the following 
discussions. 

The fiuctuation of infiaton field could be expanded into the Fourier modes as 

(^e^'-'^^0k(t) . (26) 
In this way, the fiuctuations would evolve along the equation of motion as 

<50k + 3/7(50k + A;2^/0k = , (27) 
where the effective wavenumber is given as 



e^f^ = e(i + B-t^ik-h)) , (28) 



and and k ■ n denotes the cosine of angle between the wavevector k and the third spatial 
direction fi. 

With the redefinition of the inflaton 50k = (5crk/a, we could work in the conformal time 
dr = dt/a where the coefficient (1 + i?)^ is discarded since it does not affect the following 
discussions. Thus, the equation of motion of the fluctuations becomes 

For the inflationary phase of the universe, the scale factor expands nearly exponentially, i.e., 
a{t) ~ c^*. Correspondingly, the conformal scale factor reads 

air) = --i- (r < 0) . (30) 

Thus, we flnd that the equation becomes as 

^ . U„ -^)sa.^O. (31) 



This equation has an exact solution 63] 



^ik^ffT 



5a^ = —=={l + - . (32) 



By considering the redeflnition of 50k; we obtain 

o 1 1 



1 + 



«^2|A;e//| V (l^e// 

2^2 



r 



= ^ I' + Gfc^^TFF J • ^^^^ 

The primordial power spectra of 50k, which is denoted by Vs^^, is deflned by 

27]-2 

<O|50*50^no>='5^'nk-k')— . (34) 



Thus, one has the primordial power spectra of 50k as 



f^y7Ar.f^VA(-r (35) 



\2'K J \\keff\J V27r/ \Kff\ 
For the super-horizon perturbations, the wavelength k is much larger than the horizon, i.e., 
A;r ^ 1. Thus, the primordial power spectra of 50k approximately becomes 



V27ry \\keff\ 
10 



Note that Ikeffl"^ = k'^ (^{1 + B)"^ + (3-B'/2A;)^(k ■ n)^j contains the anisotropic properties. 

On the super-horizon scales, the primordial power spectra of the comoving curvature 
perturbation 71 is given by [71 1 



where rupi is the Planck mass. The above equation could be parameterized as 

Vn{k, z) ^ r-ik) (^-W + 6B' - ^ ' (k ■ n)2 j , (38) 

where V'!^°{k) = {k / H)"''^"'^ , An denotes the normalized amplitude and uti denotes the 
spectral index of the comoving curvature perturbation. We remain only the terms that are 
not higher than i?^(r ■ fi) orders. Here n denotes the third spatial direction and r is the 
spatial location. The terms containing B induce the anisotropic behaviors of the primordial 
power spectra of the comoving curvature perturbation. The reason is that B = B{r ■ h) 
depends on the third spatial coordinate z = r ■ n. This term would give arise to the 
(parity violating) statistical anisotropy of the universe. It was noteworthy that the level 
of the statistically anisotropic effects grows with the increase of the third spatial distance 
\z\. Thus, the statistical anisotropy would be significant at large scales. At small scales, 
however, the anisotropic effects are not significant in this model. Note that the above 
primordial power spectra f l38l) could be constrained by the Planck's results via a similar 
analysis in the reference [7^,173]. 

As an example, we refer to the linear approximation of B{z) in ( 12T1) . i.e., B{z) = fiz + 
0(2;^). Here /i is a constant coefficient with the dimension [length]"^ . Thus, B becomes 

B{z)=f,r-h. (39) 



By substituting the above equation into fl38|) . one could obtain the primordial power spectra 
as 

VniK r) = ri^"{k) (^1 - 3/in ■ r + 6/1^(11 . r)^ - M(n . k)^^ . (40) 

The last two terms refer to the quadrupole asymmetry in the above equation. Generically, 
they are smaller than the dipolar term. Thus, one could discard them if one is only interested 
in the dipolar modulation of the power spectra. In this way, the primordial power spectra 
f l40|) could be rewritten as 

P^(k,r)=P-°(A;)(l-3An-r) . (41) 
11 



This dipolar modulation of primordial power spectra could be used to account for the hemi- 
spherical asymmetry (and the parity asjTiimetry) of the CMB temperature-fluctuation am- 
plitude observed by the WMAP and Planck satelhtes. 



VI. CONCLUSIONS AND DISCUSSIONS 

In this paper, we have studied a generalized FRW metric with weak anisotropy in the 
Randers spacetime. The osculating Riemannian approach was employed to obtain the evolu- 
tion of the Randers spacetime. We found an inflationary solution of Einstein's gravitational 
field equations at zero order. Meanwhile, we also obtained the anisotropic modifications to 
the inflationary phase of the very early universe. Most importantly, the primordial power 
spectra of the scalar perturbation (and the comoving curvature perturbation) was acquired 
by analysis on the equation of motion for the primordial fluctuations. It was found to be 
statistically anisotropic at large scales along the third spatial axis. This could reveal the 
deviations from the isotropy of the universe. The anisotropic effects are significant at large 
scales while they are insignificant at small scales in this model. This is consistent with the 
Plank 2013 results [2j released recently. 

The primordial power spectra with direction dependence may account for certain anoma- 



lies of the CMB observations by the WMAP and Planck satellites |74i-l80[|. In this paper, we 
obtained such a primordial power spectra with weak anisotropy in the Randers spacetime. 
Actually, the Randers spacetime has an extra vector field in the Randers structure. It could 



be related with certain (primordial) electromagnetic field with a privileged orientation [54 1. 
It singles out a preferred direction in the spacetime background. This could influence the 
inflationary evolution of the very early universe and induce the asymmetry of the primordial 
power spectra. From the primordial power spectra of the scalar perturbation in this paper, 
we could see that the parity-odd direction-dependence (such as the dipolar modulation) 
may reveal the hemispherical power asymmetry and the parity asymmetry in the CMB data 
such as the WMAP and Planck. In addition, the parity-odd part has been pointed out to 



explain the alignment of low multipoles of the CMB recently [79|, l8C [. It is noteworthy that 
the Randers spacetime could also explain the large-scale bulk flow 47|] and the anisotropic 
indications of the Hubble diagram |51| . 

As we have mentioned above, Finsler geometry is a reasonable candidate platform to 

12 



study the spacetime asymmetry and anisotropy. The Finsler spacetime is intrinsically 
anisotropic and direction dependent. The physical processes in Finsler spacetime should 
acquire similar anisotropic properties. In this paper, we indeed obtianed the anisotropic 
effects on the inflation and the primordial power spectra from Finsler geometry. To our 
knowledge, this is the first time that the primordial power spectra was obtained in the 
Finsler geometric framework. Meanwhile, the statistically anisotropic properties were pre- 
dicted for the primordial power spectra. Fortunately, the Planck's observations showed 
evidence for the deviations from isotropy. This is a good chance to test the anisotropic 
predictions of Finsler geometry. Inversely, Finsler geometry provides a natural explanation 
to the statistical anisotropy of the universe. 



ACKNOWLEDGMENTS 

We would like to thank Yunguo Jiang, Miao Li, Ming-Hua Li, Xin Li and Hai-Nan Lin 
for helpful discussions. Meanwhile, we are grateful to Banning Li for his help in our analytic 
calculations. Thanks David Lyth for his useful comment. This work has been funded by the 
National Natural Science Fund of China under Grant No. 11075166 and No. 11147176. 



[1] Planck Collaboration, larXiv: 1303.50621 

[2] Planck Collaboration, larXiv: 1303.50831 

[3] M. Tegmark, A. de Oliveira-Costa, and A. Hamilton, Phys. Rev. D 68, 123523 (2003). 

[4] H. K. Eriksen, F. K. Hansen, A. J. Banday, K. M. Gorski, and P. B. Lilje, Astrophys. J. 605, 
14 (2004). 

[5] F. K.Hansen, A. J. Banday, and K. M. Gorski, MNRAS 354, 641 (2004). 

[6] K. Land and J. Magueijo, Phys. Rev. D, 72, 101302 (2005). 

[7] J. Kim and P. Naselsky, Astrophys. J. 714, L265 (2010). 

[8] P. Naselsky, W. Zhao, J. Kim, and S. Chen, Astrophys. J. 749, 31 (2012). 

[9] J. Kim and P. Naselsky, Phys. Rev. D 82, 063002 (2010). 

[10] A. Gruppuso, et al., MNRAS 411, 1445 (2011). 

13 



[11] M.J. Longo, Phys. Lett. B 699, 224 (2011). 

[12] P.K. Aluri, and P. Jain, MNRAS 419, 3378 (2012). 

[13] E. Komatsu et al. [WMAP collaboration], Astrophys. J. Suppl. 192, 18 (2011). 

[14] C. Dvorkin, H. Peiris andW. Hu, Phys. Rev. D 77, 063008 (2008). 

[15] N.E. Groeneboom and H.K. Eriksen, Astrophys. J. 690, 1807 (2009). 

[16] N.E. Groeneboom, L. Ackerman, I.K. Wehus, and H.K. Eriksen, Astrophys. J. 722, 452 (2010) 

[17] C.G. Boehmer and D.F. Mota, Phys. Lett. B 663, 168 (2008). 

[18] L. Ackerman, S. Carroll, and M. Wise, Phys. Rev. D 75, 083502 (2007). 

[19] B. Himmetoglu, C.R. Contaldi, and M. Peloso, Phys. Rev. D 79, 063517 (2009). 

[20] B. Himmetoglu, C.R. Contaldi, and M. Peloso, Phys. Rev. Lett. 102, 111301 (2009). 

[21] S. Yokoyama and J. Soda, JCAP 0808, 005 (2008). 

[22] M. Watanabe, S. Kanno, and J. Soda, Phys. Rev. Lett. 102, 191302 (2009). 

[23] S. Kanno, J. Soda, and M. Watanabe, JCAP 0912, 009 (2009). 

[24] M. Watanabe, S. Kanno, and J. Soda, Prog. Theor. Phys. 123, 1041 (2010). 

[25] M. Watanabe, S. Kanno, and J. Soda, Mon. Not. Roy. Astron. Soc. 412, L83 (2011). 

[26] C.-H. Wang, Y.-H. Wu, and S.D.H. Hsu, Phys. Lett. B 713, 6 (2012). 

[27] D.H. Lyth and M. Karciauskas, l arXiv: 1302.73041 

[28] P. Finsler, Uher Kurven and Flachen in allgemeinen Rdumen, (Dissertation, Gottingan, 1918), 

Birkhauser Verlag, Basel, 1951. 

[29] H. Rund, The Differential Geometry of Finsler Spaces, Springer, Berlin, 1959. 

[30] D. Bao, S.S. Chern, and Z. Shen, An Introduction to Riemann-Finsler Geometry, Graduate 

Texts in Mathmatics 200, Springer, New York, 2000. 

[31] Z. Shen, Lectures on Finsler Geometry, World Scientific, Singapore, 2001. 

[32] H.C. Wang, J. London Math. Soc. sl-22 (1), 5 (1947). 

[33] S.F. Rutz, Contemp. Math. 196, 289 (1996). 

[34] X. Li and Z. Chang, Differ. Geom. Appl. 30, 737 (2012). 

[35] G.Yu. Bogoslovsky, II Nuovo Cimento B 40, 99 (1977). 

[36] G.Yu. Bogoslovsky, II Nuovo Cimento B 40, 116 (1977). 

[37] G.Yu. Bogoslovsky, II Nuovo Cimento B 43, 377 (1978). 

[38] J.R. Ellis, N.E. Mavromatos and D.V. Nanopoulos, Phys. Rev. D 61 (1999) 027503. 

[39] F. Girelh, S. Liberati and L. Sindoni, Phys. Rev. D 75 (2007) 064015. 

14 



[40] G.W. Gibbons, J. Gomis and C.N. Pope, Phys. Rev. D 76 (2007) 081701. 

[41] V.A. Kostelecky, Phys. Lett. B 701 (2011) 137. 

[42] Z. Chang and S. Wang, Eur. Phys. J. C 72, 2165 (2012). 

[43] Z. Chang and S. Wang, Eur. Phys. J. C 73, 2337 (2013). 

[44] A.P. Kouretsis, M. Stathakopoulos, and P.C. Stavrinos, Phys. Rev. D 79, 104011 (2009). 
[45] A. Kashhnsky, F. Atrio-Barandela, D. Kocevski, and H. Ebehng, Astrophys. J. 686, L49 
(2009). 

[46] A. Kashhnsky, F. Atrio-Barandela, H. Ebehng, A. Edge, and D. Kocevski, Astrophys. J. 712, 
L81 (2010). 

[47] Z. Chang, M.-H. Li, and S. Wang, larXiv:1303.1596l 

[48] L Antoniou and L. Perivolaropoulos, JCAP 1012, 012 (2010). 

[49] R.-G. Cai and Z.-L. Tuo, JCAP 1202, 004 (2012). 

[50] B. Kalus, D.J. Schwarz, M. Seikel, and A. Wiegand, arXiv : 1212.36911 

[51] Z. Chang, M.-H. Li, X. Li, and S. Wang, larXiv: 1303. 15931 

[52] Z. Chang, S. Wang, and X. Li, Eur. Phys. J. C 72, 1838 (2012). 

[53] X. Li, M.-H. Li, H.-N. Lin, and Z. Chang, MNRAS, 428, 2939 (2013). 

[54] G. Randers, Phys. Rev. 59, 195 (1941). 

[55] C.W. Misner, K.S. Thorne, and J. A. Wheeler, Gravitation, Freeman, San Francisco (1973). 
[56] G.S. Asanov, Finsler Geometry, Relativity and Gauge Theories, Reidel, Dordrecht (1995). 
[57] P.C. Stavrinos, A.P. Kouretsis, and M. Stathakopoulos, Gen. Relativ. Gravit. 40, 1403 (2008). 
[58] A. A. Starobinsky, Phys. Lett. B 91, 99 (1980). 
[59] A.H. Guth, Phys. Rev. D 23, 347 (1981). 
[60] A.D. Linde, Phys. Lett. B 108, 389 (1982). 

[61] A. Albrecht and P.J. Steinhardt, Phys. Rev. Lett. 48, 1220 (1982). 
[62] A.D. Linde, Phys. Lett. B 129, 177 (1983). 

[63] S. Dodelson, Modern Cosmology, Elsevier Pte Ltd., Singapore, 2008. 

[64] V.F. Mukhanov, and G. Chibisov, JETP Lett. 33, 532 (1981). 

[65] V.F. Mukhanov, and G. Chibisov, Sov. Phys. JETP 56, 258 (1982). 

[66] S. Hawking, Phys. Lett. B 115, 295 (1982). 

[67] A.H. Guth, and S. Pi, Phys. Rev. Lett. 49, 1110 (1982). 

[68] A.A. Starobinsky, Phys. Lett. B 117, 175 (1982). 

15 



J.M. Bardeen, P.J. Steinhardt, and M.S. Turner, Phys. Rev. D 28, 679 (1983). 
V.F. Mukhanov, JETP Lett. 41, 493 (1985). 

A. Riotto, Lectures given at the Summer School on Astroparticle Physics and Cosmology 



(2002), arXiv:hep-ph/0210162 



Y.-Z. Ma, G. Efstathiou, and A. Challinor, Phys. Rev. D 83, 083005 (2011). 
A.R. Pullen and M. Kamionkowski, Phys. Rev. D 76, 103529 (2007). 
R.Durrer, T. Kahniashvih, and A. Yates, Phys. Rev. D 58, 123004 (1998). 
J. Kim and P. Naselsky, JCAP 7, 41 (2009). 

L. Campanelh, P. Cea, and L. Tedesco, Phys. Rev. Lett. 97, 131302 (2006). 

L. Campanelh, P. Cea, and L. Tedesco, Phys. Rev. D 76, 063007 (2007). 

S. Baghram, M.H. Namjoo, and H. Firouziahi. larXiv:1303.4368l 

C. Gordon, W. Hu, D. Huterer, and T. Crawford, Phys. Rev. D 72, 103002 (2005). 

P.K. Rath, T. Mudholkar, P. Jain, P.K. Aluri, and S. Panda, JCAP 04, 008 (2013). 



16 



